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Abstract. Many quantum control tasks aim at manipulating the state of a quantum 
mechanical system within a finite subspace of states. However, couplings to the outside 
are often inevitable. Here we discuss strategies which keep the system in the controlled 
subspace by applying strong interactions onto the outside. This is done by drawing 
analogies to simple toy models and to the quantum Zeno effect. Special attention is 
paid to the constructive use of dissipation in the protection of subspaces. 



1. Introduction 

Quantum control techniques like Hamiltonian engineering are very successful when 
manipulating finite-dimensional Hilbert spaces [H El L3]. However, quantum control 
tasks often require the control of an infinite-dimensional state space. This applies for 
example, when the quantum system of interest couples to an infinitely large reservoir 
or when it contains bosonic modes with infinitely many states. In such a situation one 
should take advantage of mechanisms which restrict the time evolution of the system 
effectively onto a finite-dimensional space. The corresponding evolution can then be 
used to engineer appropriate control sequences jU [5] . 




Figure 1. Illustration of the control problem. We wish to control the system evolution 
within a subspace (yellow) of the total space (white) whilst there is coupling to an 
external subspace (blue). 
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The protection of finite-dimensional subspaces against the leakage of population 
(cf. Fig. [1]) can be achieved in many different ways. One approach is to use Hamiltonians 
which act only on small subsets of states and naturally restrict the time evolution of the 
system onto a finite-dimensional state space. This technique is used for example in ion 
trap quantum computing, where gate operations are realised by applying the required 
interactions in successive steps [6J. During each time step, the number of phonons in 
the system increases at most by one. The populated state space remains finite and 
the excitation of coherent phonon states has been avoided. At the end of every control 
sequence, the ions return into a state with zero phonon excitation [7J [8] . 

Another very efficient tool for the protection of subspaces against leakage errors are 
stimulated Raman adiabatic passages (STIRAP) [91 [10]. These employ the adiabatic 
theorem to induce transitions between states with no direct coupling between them. In 
composite quantum systems, like atoms which move slowly through an optical cavity, 
STIRAP can create ground state entanglement without populating excited electronic 
states and without creating photons inside the resonator [UJ. An alternative approach 
for protecting subspaces against leakage errors is to use especially designed impulsive 
pulse sequences. This technique has initially been developed to minimise radiation 
damage when exciting specific ground-state vibrational modes of molecules [12J. A 
generalisation of both strategies is to simply use numerical simulations which impose 
state dependent constraints to design optimal control sequences [13] . 

In this paper we discuss a very different strategy for protecting the controlled 
subspace against leakage. Instead of imposing a well designed dynamics onto its states, 
we consider strong interactions which act only on the outside. If these interactions 
introduce a time scale into the system which is much shorter than the time scale 
on which the leakage of population out of the controlled subspace would occur in 
the unprotected case, one can show that these unwanted transitions become strongly 
inhibited [Ml [151 HB]- This approach has many similarities with bang-bang and its 
generalisation dynamical decoupling [TTJ [T8j [19] which also interrupt a relatively slow 
evolution with strong interactions. 

The second half of the paper pays special attention to the constructive use of 
dissipation in the protection of subspaces against population leakage and assumes 
non-zero spontaneous decay rates of the outside states. We are interested in cases, 
where the time evolution of the system can be understood in terms of rapidly- 
repeated measurements whether the system remains in the controlled subspace or not 
[TT1 120~1 [2TJ [22l [231 1211 [25] . These force the system to remain there much longer than 
in the unprotected case which can be understood in terms of the quantum Zeno effect 
[261 1271 [28] . Notice that this approach provides a build-in error detection mechanism, 
when it is possible to register and act upon unwanted measurement outcomes. Moreover, 
as we shall see below, dissipation can protect subspaces against leakage errors even in 
situations, where other methods would simply fail. 

Refs. [291 EH] discuss similarities between bang-bang, dynamical decoupling, the 
protection of subspaces with strong interactions, and the protection of subspaces using 
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Figure 2. Toy model illustrating the leakage of population from an unprotected 
controlled subspace (represented by |0)) with coupling strength £ into an outside space 
(represented by |1)). 

dissipation. The authors conclude that all these approaches are essentially equivalent, 
since all of them can be understood in terms of the quantum Zeno effect [26]. The 
purpose of the present manuscript is to give more insight into the underlying processes. 
This is done by analysing relatively simple toy models which allow us to compare the 
above mentioned methods qualitatively as well as quantitatively. 

This paper is organized in five sections. In the next section, we consider a two- 
level system with resonant coherent coupling to obtain information about the expected 
leakage rates in a simple unprotected subspace scenario. In Section [3J we extend the 
outside space and show how strong interactions acting on the outside space can be used 
(or how they should not be used) to protect the controlled subspace against leakage 
error. Section H] analyses closely related level schemes but with non-zero spontaneous 
decay rates. An example is given where dissipation results in the protection of a subspace 
which would not be there otherwise. We finally summarise our findings in Section [5j 

2. An unprotected subspace 

Let us first consider a case where no effort is made to protect a controlled subspace from 
leaking population into outside states. For simplicity, we assume that the controlled 
subspace contains only a single state |0). As shown in Fig. [2} there is moreover only 
one relevant state outside the controlled subspace which we denote by The leakage 
of population from level into level 1 could be due to resonant interactions (like a laser 
field). Although this is an almost trivial case, the analysis of the time evolution of this 
level scheme introduces the relevant time scales of the system. This will enable us later 
to characterise and to compare the effectiveness of different strategies for the protection 
of controlled subspaces against leakage errors. 

In the following, we assume that the laser is in resonance with the 0-1 transition 
and denote its (real) Rabi frequency by £. Moreover, hcoi denotes the energy of states 

Then the system Hamiltonian in the usual rotating wave and dipole approximation 
can be written as 

i 

H = h^e^ 1 0> < 1 1 + h.c. +^2hUi (1) 

8=0 

with uj£ = (jJi — ujq. To solve the corresponding time evolution, we first change into an the 
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(5) 



interaction picture with respect to H = Yli=o ^* This transfers the Hamiltonian 

([1]) into the interaction Hamiltonian 

#i = ££|0)(l|+h.c. (2) 

To estimate the leakage rate of the controlled subspace in this case, we now calculate 
the population P m |0) at time t, given that the system was initially in |0). 
One way of doing this is to consider the usual Pauli operators o~2 and 03, 

a 2 = -i(|0)(l|-|l)(0|) and a 3 = |0)(0| - |1)(1| , (3) 

and to use the relation 

(A)= -^([AH,]) (4) 

for the time evolution of the expectation value of an operator A in the interaction picture 
to obtain a closed set of rate equations. This yields the differential equations 

which can be solved easily analytically. 

A more straightforward way of solving the time evolution of the system is to write 
its state vector as \if)) = J2 i=0 1 Ci \i) and to use the Schrodinger equation to obtain 
differential equations for the complex coefficients Cj. However, the above approach 
of deriving rate equations for expectation values is more efficient, since we are only 
interested in the leakage of population out of the controlled subspace. Moreover, this 
approach can be extended easily to include more complex level schemes as well as the 
effect of spontaneous photon emission. 

Since 0-3 commutes with H Q , we can calculate Pq using the relation 

P (t) = \(l + Mt)}) ■ (6) 

Solving Eq. (J5J) for time-independent coupling constants £ and for the case where the 
system is initially in |0), we find that the population in the initial state changes according 
to 

P Q (t) = i (1 + cos(2£t)) = cos 2 (£t) . (7) 

This means, in the absence of any protection, the system remains inside the controlled 
subspace only on a time scale which is short compared to l/£. 



3. Protecting a subspace with strong interactions 

One way to protect the controlled subspace against errors is to involve the relevant 
outside states into a relatively fast time evolution. Indeed it has been found that 
strong interactions can have the same effect as rapidly repeated measurements whether 
the system remains in its initial subspace or not [2§J (3D]- In good agreement with 
the predictions of the quantum Zeno effect [26J, these measurements strongly inhibit 
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Figure 3. Toy models to illustrate the possible protection of a controlled subspace 
(represented by |0)) with strong interactions with coupling strength f2 in the outside 
space. Here the outside space contains either the two states |1) and |2) (a) or the three 
states |1), 1 2), and |3) (b). 

transitions out of the controlled subspace. In the following, we illustrate this approach 
with the help of the two toy models shown in Fig. |3j The purpose of the interactions 
with amplitude Q is to induce fast oscillations of the amplitude of the state |1). These 
cause (a 2 ) in Eq. (J3J) to oscillate rapidly in time, such that (d 3 ) in Eq. (jSJ) becomes zero 
on average and the system remains approximately in |0). 

As we shall see below, this strategy works well for the level scheme in Fig. [3^ a). 
However, strong interactions acting on the outside space do not always protect the 
controlled subspace against leakage errors. Problems arise for example in the level 
scheme in Fig. [3](b) . There the interactions in the outside space are more complex 
than the interactions which cause the leakage. The result is that the generation of 
approximate dark states in the outside space. These are zero eigenvectors of the fast 
system dynamics. Transitions between dark states and the controlled subspace are hence 
not protected by time scale separation, even when Q becomes very large. 

3.1. Single- coupling case 

We begin with an analysis of the three-level system shown in Fig. [3]^a). Again, the 
controlled subspace contains only a single state, |0), while the outside subspace contains 
the two states |1) and |2). In order to maximise the effect of the applied interactions, we 
assume resonant couplings. As before, £ is the coupling constant for the 0-1 transition, 
while Vt denotes the coupling constant for the 1-2 transition. Here we are especially 
interested in the case, where £ <C Q. Again we have a closer look at the time evolution 
of the population P in the controlled subspace. 

As in Section |2j we denote the energy of level i by hoji. The Hamiltonian for the 
level configuration in Fig. E^a) can then be written as 

2 

H = h^e^ 1 0> < 1 1 + hrie i0Jnt |1)(2| + h.c. (8) 

j=0 
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with uj^ = oji — loq and ujq = U2 — U\. Transforming this Hamiltonian into the interaction 

Y^i=o we obtain 

-hn\i)(2\ + h.c. 



picture with respect to if 
#i = ^|0)(l 



(9) 



This interaction Hamiltonian is time independent and contains only the weak coupling 
between |0) and |1) and the strong coupling between |1) and |2). 

In order to obtain a closed system of rate equations, we now consider the expectation 
values of the Gell-Mann matrices ETJ 



a 2 
0-3 



-i(|0)(l 
[0><0| — 1 



1)(0|) • 

0"8 



0-7 



= -i(|l)(2|-|2)(l|), a 4 
(|0)(0| + |1)(1|-2|2)(2|) 



|0)(2| + |2)(0| 



(10) 



These are generalisations of the Pauli operators used in Section [2j Overall there are eight 
Gell-Mann matrices which can be used to model the time evolution of coupled three- 
level systems in a convenient way However, due to the specific form of the interactions 
in the level scheme in Fig. E(a), we need to consider only five of them. Using relation 
(j3J), we find the following closed system of differential equations 



These differential equations can be solved for example by calculating analytical 
expressions for the eigenvalues of this matrix or by simply using Mathematica. 
From Eq. ffTO]) we see that the population in the controlled subspace equals 
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(12) 



Substituting the solution of the above rate equations into this equation, we find that 
the population in |0) at time t evolves according to 



Po(t) 



2VL A + £ 4 2Q 2 ^ 2 



2/i 4 



+ 



£ 4 

cos(/it) + — — cos(2//t) 



/x 4 ^ ' 2/i 4 



(13) 



with fi 2 = Q 2 + £ 2 , if the system was initially in |0). For £ <C Q, Eq. ( Tl3l) simplifies to 
Po(t) = l-%[l-cos(nt)] (14) 



which holds up to first order in £ 2 /f2 2 . As shown in Fig. HI the system remains to a very 
good approximation, i.e. up to variations with an amplitude proportional to £ 2 /^ 2 , in 
|0). This means, for £ 2 <C Q 2 , the controlled subspace is effectively protected against 
leakage errors. 
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time (ji/S=) 

Figure 4. Time dependence of Pq for the level scheme shown in Fig. [3ja) for different 
ratios of The system is initially in |0). For O = 0, the system leaves its initial 

state space on a time scale given by l/£. For fl > 10 £, the system remains there with 
a fidelity above 95 % which constitutes an effective protection of the initial state space. 



3.2. Double- coupling case 

Using the same notation as in the previous subsection, the Hamiltonian for the level 
configuration in Fig. 131(b) in the Schrodinger picture equals 

3 

H = ^e iw «'|0)(l| + hQj Wnt (|1)(2| + |2)(3|) + h.c. + ^/iw< . (15) 

»=o 

Again we first simplify this Hamiltonian by changing into the interaction picture with 
respect to the free Hamiltonian H = ^2 i=0 fu^i This yields 

Hi = |0)(1| + hQ (|1)(2| + |2)(3|) + h.c. (16) 

Instead of solving the corresponding Schrodinger equation, we apply again Eq. OH) to 
obtain a closed set of rate equations. 

To predict the time evolution of the population P in the controlled subspace, we 
now have to consider nine generalised Gell-Mann matrices [31] . These are 

a 2 = -i(|0)(l|-|l)(0|), a 7 = -i(|l)(2|-|2)(l|), 
<7 10 = -i(|0><3|-|3><0|), a 14 = -i(|2)(3|-|3)(2|), 
cr 4 = |0)(2| + |2)(0|, a u = |l)(3| + |3)(l|, 



a 3 = |0)(0|-|l)(l|, a 8 



1 



1 



(|0)(0| + |1)(1|-2|2)(2|) 



^i5 = ^(|0)(0| + |l)(l| + |2)(2|-3|3)(3|). 

Moreover, we notice that the interaction Hamiltonian ffTBI) can be written as 
Hi = hX 0i + hVL (cr 6 + 0-13) 



;i7) 
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Figure 5. Time dependence of Po for the level scheme shown in Fig. (3^b) for different 
ratios of Q/£. Here the controlled subspace is no longer protected against leakage, even 
when f2 becomes as large as 100 £. The reason is the leakage of population into the 
dark state |Ao) which is illustrated in Fig.H 



with 

a 1 = |0)(l| + |l)(0|, a 6 = |l)(2| + |2)(l|, a 13 = |2)(3| + |3)(2| . (19) 



Substituting Eqs. (TT7l) -(fT9l into Eq. (@J and evaluating the relevant commutators, we 
see that the expectation of the operators in Eq. (JTTJ) evolve according to 
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This system of linear differential equations can, in principle, be solved analytically. 
However, for simplicity, we restrict ourselves to the presentation of a numerical solution. 

Using the Gell Mann matrices defined in Eq. (I17p . the population in the initial state 
Po can now be written as 

p ° = i + ^> + ^> + 5vI^>- < 21 > 

The time evolution of Po obtained from substituting the numerical solution of the 
differential equation f l20|) into Eq. (121]) is shown in Fig. |5j Comparing the result for 
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Figure 6. Illustration of the effect of the Hamiltonian (|23p onto the states involved 
in the time evolution of the system. This level scheme is identical to the one shown in 
Fig. [3th) but now we clearly see why the initial state |0) is no longer protected against 
leakage errors. 



different values of 0/£ with the time evolution in the Q = case, we see that the 
controlled subspace is not protected, even when Q is much larger than £. Leakage of 
population out of the controlled subspace happens on the same time scale as in the 
unprotected case. 

Why does the protection of the controlled subspace work in the level scheme shown 
in Fig. [3]^a) but not in the very similar level scheme shown in Fig. [3](b) ? The reason for 
this becomes clear when we rewrite the Hamiltonian in Eq. (fl6|) in terms of the states 
|0>, 

|A ) = -^(|1)-|3», |A a > ^ ^(Il> + |3» (22) 
and 1 2). Using this notation, H\ becomes 

Hj = -^hi |0)(A | + 4=^ |0)(Ai| + V2MI \X 1 )(2\ + h.c. (23) 
v 2 v 2 

The effect of this Hamiltonian is illustrated in Fig. |6j It shows that the system is 
only protected against leakage into the |Ai) state, since this state experiences a strong 
interaction. However, the system is not protected against leakage into |A ), since |A ) 
is a zero eigenstate of the Q terms in Eq. fl23|) . This means, | Aq) is not involved in a 
fast evolution and the transfer from |0) to |Ao) occurs on the same time scale as in the 
unprotected case. In the final section, we show that dissipation is able to remove such 
dark states from the system so that the controlled subspace becomes protected again. 



4. Protecting a subspace with dissipation 

In this section we analyse three examples (cf. Fig. [7j) which illustrate the possible 
protection of the controlled subspace using dissipation. The controlled subspace contains 
again only the |0) state, while the outside space contains one, two or three states. The 
only difference to the examples discussed in Sections |2] and [3] is the presence of a non- 
zero spontaneous decay rate T. As we shall see below, the controlled subspace is well 
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(a) 



(b) 
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n 



Figure 7. Toy models to illustrate the possible protection of the controlled subspace 
(represented by |0)) with a non-zero spontaneous decay rate T and strong interactions 
with coupling strength D, in the outside space. 



protected against leakage in all three scenarios, when the interactions in the outside 
space described by Q and the spontaneous decay rate T are sufficiently larger than £. 



4-1. Single- coupling case with dissipation 

Let us first have a look at the three- level system shown in Fig. [7(b). To describe its time 
evolution, we go again into the interaction picture with respect to the free evolution and 
analyse the master equation 



P = ~j:[Hi,p] + \ 2 |1)(2| p |2)(1| - p |2)(2| - |2)(2| p 



(24) 



The interaction Hamiltonian H\ is the same as in Eq. ([9]). In order to predict the time 
evolution of the population in the controlled subspace, we derive again a closed system 
of rate equations. The time derivative of the expectation value of an operator A is now 
given by 

(A) = Tr (Ap) . (25) 

Taking this into account, we find that the Gell Mann matrices in Eq. ( 1101) evolve 
according to 



(26) 



These equations resemble the ones shown in Eq. ( jTlj) . The additional T terms take the 
effect of dissipation into account. 

The population in the controlled subspace can be obtained by substituting for 
example the numerical solution of these equations into Eq. ( fT2i) . The result is shown in 
Fig. |8j For simplicity we assumed r = Q. For Q = 0, we see again the Rabi oscillations 
of the unprotected case. However, when Q becomes sufficiently larger than £, then 
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Figure 8. Time dependence of Pq for the level scheme shown in Fig. [TJb) for T = 17 
and different ratios of f2/£. On average, the protection of the controlled subspace is 
more or less the same as in Fig. [4] which corresponds to the same level scheme but with 
r = (cf. Fig. Eta)). 

the system remains to a very good approximation in its initial state. On average, the 
protection of the controlled subspace is more or less the same as in Section I3TT] where we 
had T = (cf. Fig. H]). There seems to be no advantage of having a non-zero spontaneous 
decay rate in the system! Notice that having F" 7^ in the level scheme in Fig. E^b) is 
only advantageous when someone actually observes whether the system emits photons 
or not. Indeed, one can show that the system remains in its initial state |0) with a very 
high fidelity under the condition of no photon emission [32]. If a photon emission is 
detected, then the system has left the controlled subspace and the anticipated control 
experiment needs to be restarted. 

Comparing the level scheme in Fig. UXb) with the level scheme analysed in 
Refs. [32j [331 El], we see that its dynamics exhibits so-called macroscopic light and 
dark periods. Indeed, for £ much smaller than Q and T, the initial state |0) is an 
approximate zero eigenstate of the system dynamics. The absence of photon emissions 
hence confirms that the system is in this state. As a consequence of the quantum 
Zeno effect, it therefore remains there for a relatively long time. On average, this time 
equals f2 2 /r£ 2 which is much larger than l/£ [32]. In other words, the system exhibits 
a macroscopic dark period. The system may eventually drop out of the controlled 
subspace, thereby entering a so-called macroscopic light period and causing fluorescence 
at a rate which depends on Q and T. This behaviour is not reflected in Fig. [8j since 
the density matrix description used in this paper does not allow us to distinguish the 
different trajectories of the system. 
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time 

Figure 9. Time dependence of Pq for the level scheme shown in Fig. (7{a) for different 
ratios of T /£. For r>^, the system remains in the controlled subspace with a very 
high fidelity. 



4-2. Single-state outside 

Let us now have a look at the simple level configuration in Fig. [7]^a). Its time evolution 
is given by the master equation 



~[F I ,p] + |[2|0)(l|p|l)(0| 



(27) 



(28) 



i> • ~ T . //,./;]• •- ^>MJ|/>|i,HU|-p|l)(l|-|l)(l|p 

In the interaction picture with respect to the free evolution, the interaction Hamiltonian 
Hi is the same as in Eq. ([2]). To predict the time evolution of Pq we proceed as in Section 
[2] and consider again the Pauli operators in Eq. ([3]). Their expectation values evolve 
now according to 

J " \ ^ -E ){ (as) 

Combining the result of this equation with Eq. ([6]) yields the time dependence of the 
population P in the controlled subspace. 

Fig. [9] shows a numerical solution of the time dependence of Pq for different ratios 
r/£. For T = 0, we observe the Rabi oscillations in and out of the initial subspace 
which occur in the unprotected case. For V 3> £, the state vector becomes |0) with a 
very high fidelity. But even for relatively modest values for r/£, the density matrix p 
settles quickly into a steady state with the system predominantly in |0). The reason for 
this very strong protection of the controlled subspace is that, even when it leaves, the 
system returns very rapidly via the spontaneous emission of a photon. 




4-3. Double- coupling case with dissipation 

This final subsection analyses the time evolution of the four-level system shown in 
Fig. E(c). The only difference to the level scheme in Fig. E^b) is the presence of the 
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£2/g= 1 
Q/|= 10 
Q/| = 100 
Q/E = 1000 



time (nfc) 



Figure 10. Time dependence of Pq for the level scheme shown in Fig. (7{c) for 
r = ft and different ratios of f2/£. Compared to Fig. [5J we now observe an increasing 
effectiveness of protection of the controlled subspace with increasing values of Q. 



non-zero spontaneous decay rate T. To calculate the time evolution of the population 
in the controlled subspace, we now consider the master equation 



P 



r 

+ 2 



2|l)(2|p|2)(l|-p|2)(2| 
2|2)(3|p|3)(2|-p|3)(3|-|3)(3|p 



|2)(2|p 



(29) 



161) . Proceeding as above, we find 
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Fig. [TU] shows the time dependence of Pq for the case where the system is initially in 
the controlled subspace and has been obtained by substituting the numerical solution 
of these equations into Eq. ( 1211 . 

Comparing Figs. [5] and [TQl we see that the presence of a sufficiently large 
spontaneous decay rate T combined with the presence of a relatively large coupling 
constant Q now results in an effective protection of the controlled subspace against 
leakage errors. There are different ways of seeing how this protection (which was not 
there before) has been achieved. One way is to have a closer look at the above master 
equation and to notice that the state |Ao) is no longer a zero eigenstate of the system 
dynamics. Whenever, population accumulates in this state, the system returns (either 
via the emission of a photon or as a result of its no-photon evolution) on the time scale 
given by T into |1), where it experiences fast driving with Q. This example confirms 
that dissipation can provide a very efficient tool for restricting the time evolution of a 
system onto a controlled subspace. 

Another way to gain an intuition into the behaviour of the level scheme in Fig. [7(c) 
is to compare it to the level scheme in Fig. [7J(b) which we analysed in Section 14.11 
Observing whether the system emits photons or not, one would notice again two very 
distinct phases of operation. The system either emits photons at a high rate or it remains 
dark for a relatively long time. A macroscopic light period, on one hand, indicates that 
the state vector lies entirely outside the controlled subspace. A macroscopic dark period, 
on the other hand, indicates that the system is in |0). In other words, if the system 
is initially in the controlled subspace, it remains there on average much longer than in 
the unprotected case. The result is the protection of the system against leakage errors 
which, when they occur, are heralded by an easy-to-detect fluorescence signal. 

5. Conclusions 

This paper illustrates two methods to protect a controlled subspace against the leakage 
of population into the outside space: one using strong interactions in the outside space 
and one using dissipation. This is done with the help of relatively simple toy models 
whose time evolution can be analysed relatively easily. For simplicity, we assume that 
the controlled subspace consists only of one state, namely |0). The outside space contains 
either one, two or three states denoted |1), |2), and |3) (cf. Figs. |3]and[7J- Section [2] 
discusses the unprotected case and shows that unwanted transitions from |0) to |1) (due 
to resonant coupling) occur on a time scale given by a relatively small parameter £ 

(cf. Fig. ED. 

In Section [3j the decoherence time of the system is increased to one which scales as 
£ 2 by applying relatively fast interactions with coupling strength fl to the outside space. 
However, these strong interactions are not always sufficient for protecting a controlled 
subspace against leakage errors. While it works well for the level scheme shown in 
Fig. E](a), no protection occurs for the level scheme shown in Fig. [3(b). The reason is 
the existence of an approximate zero eigenstate outside the controlled subspace. This 
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state does not experience fast driving and therefore behaves as the state |1) in the 
unprotected case. 

Section @] considers three scenarios where a spontaneous decay rate T has been 
added to the level schemes analysed in Sections [2] and [31 All the level schemes shown in 
Fig. [7] exhibit a strong protection of the controlled subspace. One way to understand the 
mechanism which inhibits the population transfer out of the controlled subspace is to 
interpret the behaviour of the system in terms of the quantum Zeno effect [26, 27, 28J. 
Suppose being outside the controlled subspace results necessarily in the spontaneous 
emission of a photon. Then, observing whether a photon emission takes place or not 
is equivalent to performing a measurement on whether the system is in the controlled 
subspace or not. If these measurements occur on a sufficiently short time scale, then 
a system initially in the controlled subspace remains there much longer than in an 
unobserved case. A similar interpretation applies to the protection of the controlled 
subspace with a strong interaction illustrated in Fig. |3](a) [291 130] . 

Finally, let us remark that the absence of decoherence within the controlled 
subspace results in general in an effective time evolution which can be described by 
the effective Hamiltonian [HI \T7\ [22| [35] 

# eff = P C s#Pcs, (32) 

where Pes denotes the projector onto the controlled subspace and H is the total system 
Hamiltonian. Once the protection is in place, the interactions described by H e g can be 
designed as required by the control task at hand. Applications of such control tasks can 
be found for example in quantum information processing, i.e. in the realisation of gate 
operations and the preparation of highly entangled states. Recent ideas for achieving 
these tasks use dissipation in an even more constructive way, for example, by heralding 
successful state preparations with macroscopic fluorescence signals [361 E3 EEJ [39] 
or by letting non-unitary evolutions guide the system into the desired target states 
dQl SH S21 SHI SH SH]. 
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